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We present recent results on noise-induced transitions in a nonlinear oscillator with randomly 
modulated frequency. The presence of stochastic perturbations drastically alters the dynamical 
behaviour of the oscillator : noise can wash out a global attractor but can also have a constructive 
role by stabilizing an unstable fixed point. The random oscillator displays a rich phenomenology but 
remains elementary enough to allow for exact calculations : this system is thus a useful paradigm 
for the study of noise-induced bifurcations and is an ideal testing ground for various mathematical 
techniques. We show that the phase is determined by the sign of the Lyapunov exponent (which 
can be calculated non-perturbatively for white noise), and we derive the full phase diagram of the 
system. We also investigate the effect of time-correlations of the noise on the phase diagram and 
show that a smooth random perturbation is less efficient than white noise. We study the critical 
- - - behaviour near the transition and explain why noise-induced transitions often exhibit intermittency 

o : and multiscaling : these effects do not depend on the amplitude of the noise but rather on its 

' power spectrum. By increasing or filtering out the low frequencies of the noise, intermittency and 

' multiscaling can be enhanced or eliminated. 
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I. INTRODUCTION 



Most patterns observed in nature are created by instabilities that occur in an uncontrolled noisy environment : 
convection in the atmospheric layers and in the mantle are subject to inhomogeneous and fluctuating heat flux: sand 
dunes are formed under winds with fluctuating directions and strengths. These fluctuations usually affect the control 
(-H , parameters driving the instabilities, such as the Rayleigh number which is proportional to the imposed temperature 
' ^ ' gradient in natural convection. These fluctuations act multiplicatively on the unstable modes. In the same spirit, the 
', evolution of global quantities, averaged under small turbulent scales, can be represented by a nonlinear equation with 
fluctuating global transport coefficients that reflect the complexity at small scales. For instance, it has been shown 
that the temporal evolution of the total heat flux in rotating convection can be described by a non-linear equation 
with a multiplicative noise p^. The dynamo instability that describes the growth of the magnetic field of the earth 
' and the stars because of the motion of conducting fluids in their cores, is usually analyzed in similar terms : the 
: magnetic field is expected to grow at large scale, forced by a turbulent flow. Here again, the parameters controlling 
' the growth rate of the field are fluctuating . 

[ Since the theoretical predictions of Stratonovich (3], and the experimental works of Kawaboto, Kabashima and 
■ Tsuchiya 0] , it is well known that the phase diagram of a system can be drastically modified by the presence of noise 
', @i 0- Because of stochastic fluctuations of the control parameter, the critical value of the threshold may change 
C . and noise can delay or favor a phase transition 0, d, [13] ■ When a physical system subject to noise undergoes 
bifurcations into states that have no deterministic counterparts, the stochastic phases generated by randomness have 
specific characteristics^ (such as their scaling behavior and the associated critical exponents) allowing to define new 
universality classes {11.]. 

^ . , A straightforward approach to study the effect of noise on a bifurcation diagram would consist in analyzing the 
' nonlinear Langevin equation that governs the system. However, the interplay between noise and nonlinearity results 
in subtle effects that make nonlinear stochastic differential equations hard to handle 0, [H, [l^ ■ One of the simplest 
systems that can be used as a paradigm for the study of noise-induced phase transitions is the random frequency 
oscillator (for a recent and detailed monograph devoted to this subject see [13]). A deterministic oscillator with 
damping evolves towards the equilibrium state of minimal energy that represents the unique asymptotic state of the 
system. However, if the frequency of the oscillator is a time-dependent variable, the behavior may change : due 
to continuous energy injection into the system through the frequency variations, the system may sustain non-zero 
oscillations even in the long time limit if the amplitude of the noise is large enough. 

In this work, we present recent results on noise-induced transitions in a nonlinear oscillator with randomly modulated 
frequency. The online of this article is as follows. 

In section [Hi we show that an explicit calculation of the Lyapunov exponent allows to determine exactly the 
phase diagram of the noisy Duffing oscillator In the case of a single- well oscillator, a large enough noise 



destabilizes the origin (which is a fixed point for the deterministic dynamics). The double- weU oscillator presents a 
richer structure : the unstable fixed point becomes a global attractor of the stochastic dynamics when the amplitude of 
the noise lies in a well-defined interval : this is a typical example of a reentrant transition (which was observed earlier 
in the more elaborate setting of spatially extended systems [131 )■ These properties are only qualitatively modified 
when the noise has temporal correlations. 

In section Hm we discuss the scaling exponents associated with the order parameter in the vicinity of the transition. 
For a deterministic transition, it is well-known that the different types of bifurcations can be classified by some 
characteristic exponents. We show, in the stochastic case, that these exponents can be modified by the noise : the 
dynamical variable has a non-trivial scaling behaviour and its time series exhibits on-off intermittency. Following 
[H, [l^ , we explain that niultiscaling and intermittency are intimately linked with each other and that their existence 
is due to the presence of low frequency components in the noise spectrum. By using a random process with a spectral 
structure richer than that of white or Ornstein-Uhlenbeck noise, we show that the deterministic exponents can be 
recovered by a suitable low-band flitering of the noise. This fact implies that the concept of critical exponents can be 
ambiguous for noise-induced bifurcations and is certainly not as useful as it is for deterministic phase transitions. 

Section |IV| presents a discussion on perturbative perturbative expansion of the nonlinear Langevin equation. 
Whereas most studies rely on Fokker-Planck type evolution equations for the Probability Distribution Function (PDF) 
in which the noise is integrated out by mapping a stochastic ordinary differential equation into a deterministic partial 
differential equation in the phase space of the system, it is also possible to perform a direct perturbative expansion 
of the nonlinear Langevin equation by using the classical Poincare-Lindstedt method [13, [Mj . This technique has the 
advantage as compared to the more traditional approaches of being applicable to a noise with an arbitrary spectrum 
(whereas Fokker-Planck equations are only valid for white- noise) . However, due to the presence of low- frequencies 
in the noise spectrum, the perturbative expansion breaks down and diverges. This divergence is at the origin of the 
multiscaling behaviour and of the anomalous scaling exponents. The last section is devoted to concluding remarks. 



II. STOCHASTIC BIFURCATIONS OF A RANDOM PARAMETRIC OSCILLATOR 



In this section, we present the phase diagram of a nonlinear oscillator with a randomly modulated frequency. 
The case when the frequency is a periodic function of time is a classical problem known as the Mathieu oscillator 
[H, m, [131 ■ The phase diagram of the (damped) Mathieu oscillator is obtained by calculating the Floquet exponents 
(defined as the characteristic growth rates of the amplitude of the system) . This phase diagram presents an alternance 
of stable regions in which the system evolves towards its equilibrium state and of unstable regions in which, because of 
parametric resonance, the amplitude of the oscillations grows without bound. A nonlinear term is needed to saturate 
these oscillations. When the frequency of the pendulum is a random process, the role of the Floquet exponents is 
taken over by the Lyapunov exponents. The system undergoes a bifurcation when the largest Lyapunov exponent, 
defined as the growth rate of the logarithm of the energy, changes its sign. Thus, the Lyapunov exponent vanishes 
on the critical surface that separates the phases in the parameter space. This criterion involving the sign of the 
Lyapunov exponent resolves the ambiguities that were found in studies of the stability of higher moments jl^, [H, [13| 
and has a firm mathematical basis ^13j . In recent works [isl . [l6j , we have obtained the exact phase diagram of the 
random oscillator when the frequency is a Gaussian white noise. In the last part of this section, we discuss the effect 
of non- vanishing time-correlation in the noise f^sj . 



A. Instability induced by noise: the single-well oscillator 

The equation for the amplitude a; of a random Duffing oscillator with a fluctuating frequency is 

x + -fx+{uj^ + ^{t))x + Xx^ =0, (1) 

where 7 is the (positive) friction coefficient, lu the mean frequency andA the coefficient of the cubic non- linear term. 
The random fluctuations of the frequency are represented by a Gaussian white noise ^(t) of zero mean- value and of 
amplitude V 

m) = 0, 

imm) = v6{t-t'). (2) 

In this work, all stochastic differential equations are interpreted according to Stratonovich calculus. By rewriting time 
and amplitude in dimensionless units, t := Lot and a; := X^^^uj~^x, Eq. ([1]) becomes 
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FIG. 1: Phase plane plots of a typical trajectory of the noisy nonlinear oscillator. Eq. (|3]) is integrated numerically with a time 
step St = 5 10"". (a) a = 0.5, A = 20, t < 100; (b) a = 2.0, A = 20, t < 100; (c) a = 2.0, A = 60, t < 100. Early-time 
relaxation towards noisy oscillations has been omitted for clarity in cases (a) and (c). The trajectory spirals down towards the 
origin in case (b). 



where ^{t) is a delta-correlated Gaussian variable with vanishing mean value and with correlations given by 

(S(i)S(t')) = A5it~t'). (4) 

The parameters 

a = j/lu and A = V/uj^ (5) 

correspond to dimensionless dissipation rate and to noise strength, respectively. 

The equation ^ has the origin of the phase space, i; = a; = 0, as a fixed point. When the system is deterministic, 
i.e., when V = 0, the origin is a global attractor for any 7 > : for any initial condition we know that x{t) — > and 
i{t) — > when t 00. 

However, in presence of noise, the asymptotic behaviour of the system becomes more complex. In Fig. [1] we 
present, for three values of the parameters, a trajectory in the phase plane (x, v = x) characteristic of the oscillator's 
behaviour, using the numerical one step collocation method advocated in [29]. Initial conditions are chosen far from 
the origin, with an amplitude of order 1. Noisy oscillations are observed for small values of the damping parameter a 
[See Fig. [TJa]. For larger values of a the origin becomes a global attractor for the dynamics [See Fig.[TJb]. Increasing 
the noise amplitude A at constant a makes the origin unstable again [See Fig. [T]c]. 

In other words, for the stochastic oscillator, the following property is true : For a given 7 there is a critical noise 



amplitude I'c(7) such that 



V < "Dcij) x,x —^ when t oo 

DyDcil) Oscillating Stationary State 

The main problem is to calculate the value of I?c as a function of 7 and the other parameters of the problem. 

For random dynamical systems, it was recognized early on that various 'naive' stability criteria [25l.[2^. obtained by 
linearizing the dynamical equation around the origin, lead to ambiguous results. This feature is in contrast with the 
deterministic case for which the bifurcation threshold is obtained without ambiguity by studying the eigenstates of 
the linearized equations [30| . For the random oscillator, trying to determine the critical noise amplitude by studying 
the stability of the moments of the amplitude of the linearized equations leads to inconsistent results. More precisely, 
consider the random harmonic oscillator, obtained by linearizing equation ([T]) near the origin : 

x + -fx+{u}^ +^{t))x = 0. (6) 

The second moment (a:^) of the amplitude of the linearized equation converges to in the long time limit if 

<2-. 



The fourth moment {x"^) converges to if 



More generally, the moment {x'^p) is stable if 
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In fact, it was conjectured in [26| and proved in [27| that, in a linear oscillator with arbitrarily small parametric 
noise, all moments beyond a certain order grow exponentially in the long time limit. In figure [H we plot the stability 
diagram for the moments of order 2,4 and 6. We observe that the stability range becomes narrower when the order of 
the moment becomes higher. This means that any criterion based on finite mean displacement, momentum, or energy 
of the linearized dynamical equation is not adequate to show the stability of the initial problem and the bifurcation 
threshold of a nonlinear random dynamical system cannot be determined simply from the moments of the linearized 
system. In practice, the phase diagram of the non-linear random oscillator was usually determined in a perturbative 
manner by using weak noise expansions [13, HH, HU, [H, [11] . 

The reason for this failure is that in the linearized problem the statistics of the variable x is dominated by very 
large and very rare events that induce the divergence of moments of high order. These rare events are suppressed 
by the nonlinearity and the bifurcation threshold becomes well defined. We face the following paradoxes : (i) the 
stability region of the origin x = x = for the non-linear equation has a well defined threshold for any value of the 
noise amplitude : below this threshold the origin is a global attractor of the dynamics. On the contrary the linearized 
equation does not have a clear-cut bifurcation threshold, (ii) Close to the origin the non-linear term is irrelevant but 
the nonlinearity has to be taken into account to suppress the rare events that spoil the statistical behaviour of the 
system. 

The problem is to find the correct stability criterion for the non-linear equation that can be formulated as in the 
deterministic case on the linear equation and that would allow explicit calculations. 

Such a criterion does exist and is based on the Lyapunov exponent that measures the growth rate of the random 
harmonic oscillator's energy 

A = hm -(logii;) with E=^+ Lo^^ , (7) 

t^oo Zt Z Z 

where x{t) is a solution of equation ([6|). It has been shown (see [ll] and references therein) that when the Lyapunov 
exponent is negative the Fokker-Planck equation has a unique stationary solution which is the Dirac delta function 
at the origin of the phase space. This means that the origin is a stable global absorbing state. However, when 
the Lyapunov exponent is positive and extended, stationary probability distribution function exists and describes an 
oscillatory asymptotic state of the nonlinear random oscillator. These features are reminiscent of Anderson localization 
and there exists indeed a mapping between the random oscillator and a Id localization model [s^ [ssj. 




The equation of the transition line that determines the stabiHty of the origin for the nonhnear osciUator with a 
random phase is therefore given by 



A(7,c^,P) = 0. 



(8) 



The foUowing exact closed formula for the Lyapunov exponent of the system can be derived when the random frequency 
modulation is a Gaussian white noise llSl l36l 13711 : 
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(9) 
(10) 



Inserting this formula in the stability criterion ^ and using dimensionless variables we obtain the equation of the 
transition line in the (a, A) plane. For a given value of A, the critical value of the damping a = ac(A) for which the 
Lyapunov exponent vanishes is given by 

/;°°duvi^exp{-i((i-4)^+f;)} 

"c = A 7^ r/^- (11) 

The critical curve a = etc (A) is represented in Fig. [3l It separates two regions in parameter space: for a < etc (resp. 
a > ttc) the Lyapunov exponent is positive (resp. negative), the stationary PDF is an extended function (resp. a 
delta distribution) of the energy and the origin is unstable (resp. stable). For small values of the noise amplitude, 
it can be shown that the exact formula (jlip reduces at first order to the linear relation ac = ^, in agreement with 
previous perturbative calculations. For large values of the noise, we obtain ac — 0.656A^/'^ . We emphasize that the 
phase diagram drawn in Fig. [3] is exact for all values of the noise amplitude A and the damping parameter a. 



B. Stabilisation by noise: the double-well oscillator 



In a classical calculation, Kapitza (1951) showed that the unstable upright position of an inverted pendulum is 
stabilised if its suspension axis undergoes sinusoidal vibrations of high enough frequency. Analytical derivations of 
the stability limit are based on perturbative approaches, i.e., in the limit of small forcing amplitudes [H, [H,[3^. 




FIG. 3: Phase diagram of the single-weU oscillator with multiplicative noise. The critical curve separates an absorbing and an 
oscillatory phase in the (a, A) plane. The dashed line represents the small noise approximation (see text). 



When the sinusoidal vibrations of the suspension axis are replaced by a white noise, an exact, non-perturbative, 
stability analysis can be performed p^ . In particular, for the stochastic inverted pendulum, we have shown that 
the unstable fixed point can be stabilized by noise and have discovered the existence of a noise-induced reentrant 
transition. 

The general equation for the Duffing oscillator subject to multiplicative noise is given by : 



0. 



(12) 



When /.t < 0, the non-linear potential has a single-well. This case is the same as the one discussed in the previous 
section : the origin is deterministically a global attractor that can become unstable in presence of noise. When /x > 0, 
the oscillator is subject to a double-well potential and the origin is deterministically unstable. 




FIG. 4: The double-well oscillator with parametric noise 

In presence of noise, the correct criterion for stability analysis is again based on the sign of the Lyapunov exponent. 
In order to take into account both possible signs of ^, it is useful to define the following set of dimensionless parameters : 

a = and A = ^ . (13) 

An exact calculation of the Lyapunov exponent [l6| then allows to draw the phase diagram represented in figure 
We observe that when a < 0.21 . . ., there exits a range of noise amplitudes for which the origin is stabilized by 

noise : for A such that Ai(a) < A < A2(a) the origin becomes an attractive fixed point of the stochastic dynamics. 

The values Ai(a) and A2(a) that determine the stability interval are known analyticaly. We emphasize that these 
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FIG. 5: Phase diagram of the double-weU osciUator subject to muhiphcative Gaussian white noise. The solid hne is the locus 
where A(Qf, A) = 0. The bifurcation line a = of the noiseless dynamical system is drawn for comparison (dotted line). For 
A < A* ~ 3.55 (resp. A > A*), the origin is stabilised (resp. destabilised) by the stochastic forcing in the range < a < Qc(A) 
(resp. tic (A) < a < 0). 



functions can not be calculated perturbatively by using a small noise expansion when a is finite. For a > 0.21 . . ., 
the stability interval does not exist anymore : the origin is always unstable and the non-equilibrium stationary state 
exhibits an oscillatory behaviour. 
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FIG. 6: Stability range for the double- well oscillator 
noise amplitudes between Ai and A2. 



when Q < 0.21, the unstable equilibrium point is stabilized by noise for 



The effect of colored noise 



We now discuss the phase diagram of an oscillator whose frequency is a random process with finite time memory. 
More precisely, we consider the case of an Ornstein-Uhlenbeck noise x(t) of correlation time r : 



(mm) 
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exp {-\t-t'\/T) 



(14) 



When r ^ 0, the process £,{t) becomes identical to the white noise. The Lyapunov exponent of the system A(w, 7, V, r) 
becomes now a function of t also. 

From a physical point of view, the influence of a finite correlation time on the shape of the critical curve is an 
interesting open question: due to the finite correlation time of the noise, the random oscillator is a non-Markovian 
random process and there exists no closed Fokker-Planck equation that describes the dynamics of the Probability 
Distribution Function (P.D.F.) in the phase space. This non-Markovian feature hinders an exact solution in contrast 
with the white noise case where a closed formula for the Lyapunov exponent was found. For the single-well stochastic 
oscillator, we have calculated this Lyapunov exponent by using different approximations and have derived the 
phase diagram. The main role of the correlation time of the noise, as can be seen in Figure [71 is to enhance the 
stability region. When r grows, amplitude of the noise required to destabilize the origin becomes bigger. This effect 



can be seen quantitavely in the following exact asymptotic expansions (see also fioj ) 



For small values of the noise amplitude : etc 



A 



4(1 + 4^2^2) ' 



/ A \ 1/4 

For large values of the noise amplitude : Uc — C — 

\Ll!T/ 



(15) 
(16) 



where the constant C is of order 1. Comparing with the white noise result we observe that the behaviour of the critical 
curve for large noise is modified in presence of time correlations; the asymprtotic exponent is 1/3 for white noise and 
1/4 for colored noise even if the correlation time r is small. The presence of a non- vanishing correlation time thus 
modifies the scaling characteristics of the system. This change of scaling at large noise is related to the differentiability 
properties of the random potential as shown in (4]| : the white noise is continuous but nowhere differentiable whereas 
the Ornstein-Uhlenbeck process has a first derivative but no second derivative. 
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FIG. 7: Critical curves obtained by numerical simulations of the single-well random oscillator subject to multiplicative Ornstein- 
Uhlenbeck noise. The curves with black dots (•) and with squares (□) correspond to r = 1 and t — 2, respectively. The full 
black line is the analytical result for the white noise case. 

For the double-well oscillator subject to Ornstein-Uhlenbeck noise, the results of numerical computations of the 
phase diagram are shown in figure [5] Here also we observe that the noise-mducedstabilization of the deterministically 
unstable fixed point becomes less efficient in presence of time correlations. In the weak noise limit, this curve agrees 
with the prediction of ^]:ac(A) ~ A/(2(1-|-t)). For a noise amplitude of order 1, the bifurcation line is qualitatively 
similar to that obtained with white noise. However, depending on the value of A, the value of the bifurcation point 
Q:c(A) is not necessarily a monotonic function of r. A precise understanding of this non- monotonous behaviour is 
lacking and the analytical theory of the stability of the double- well oscillator with time-correlated noise still remains to 
be done. Although the obtention of exact results for the Ornstein-Uhlenbeck noise seems unlikely, various recent works 
[43 | indicate that the Poisson process provides a useful model for the study of time correlation effects in stochastic 
systems. The advantage of the Poisson process is that it is amenable to exact analytic calculations. 

We conclude this section by emphasizing that the results obtained here for the oscillator with multiplicative noise 
could be used for other stochastic systems. For example, it has been shown by Schimansky-Geier et al. [i^ that 
the random Duffing oscillator with additive noise undergoes a phase transition that does not manifest itself in the 
stationary P.D.F. (which is simply given by the Gibbs-Boltzmann formula). This subtle phase transition, which affects 
the properties of the random attractor of motion in phase space, can be formulated mathematically as a bifurcation 
in an associated linear oscillator subject to a multiplicative noise with a a finite correlation time. If we approximate 
this noise by an effective Ornstein-Uhlenbeck process, the system becomes identical to the one studied here. 



III. SCALING BEHAVIOUR NEAR THE BIFURCATION THRESHOLD 



For deterministic Hopf bifurcation, the amplitude of the order parameter x exhibits a normal scaling behaviour in 
the vicinity of the transition line; if e denotes the distance from threshold of the control parameter, we have in the 
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FIG. 8: Phase diagram of the double- well oscillator subject to Ornstein-Uhlenbeck noise with correlation time r. The bifurcation 
lines for r = 0.1 and 1.0 are obtained from numerical calculations of the Lyapunov exponent. For comparison, we draw the 
(analytic) white noise line. In this figure the dimensionless parameters are given by a = ^ and A : 



deterministic case 
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(17) 



We now consider the case of the stochastic oscillator. In figure [51 we plot the behaviour of the amphtude moments 
in the vicinity of the reentrant transition in a double-well oscillator. The damping rate a is fixed at a given value 
and the noise strength A is chosen such that < A — Ac(a) ^ 1, where Ac(a) is the critical value. In this case, 
the parameters of the system are tuned just above the bifurcation threshold and the Lyapunov exponent is slightly 
positive: < A <C 1. We observe that even-order moments of the amplitude scale linearly with the distance to 
threshold in the vicinity of the bifurcation line (note that odd-order moments are equal to zero by symmetry) : 



(18) 



A similar scaling behaviour is obtained by simulating the single- well oscillator or by replacing white noise by an 
Ornstein-Uhlenbeck process. Such a behaviour was also observed in random maps [44| | and in stochastic fields [45|. 
Such a strong multiscaling seems to be a typical feature of stochastic bifurcations and has been noticed a long ago [5| 
in first order stochastic differential equations. The fact that the stochastic oscillator exhibits a similar behaviour can 
be understood 'physically' by noticing that the second order time derivative becomes irrelevant in the long-time limit. 
The mathematical explanation for this behaviour is also elementary : near the bifurcation threshold, when the 
Lyapunov exponent A — » 0, the stationary Probability Distribution Function Pgtat exhibits a power law divergence at 
the origin {x = x = 0), and the region near the origin dominates the statistics : it can be shown that near E — 0, 
the stationary energy distribution scales as Pstat{E) ~ E'^^~^ where c is a positive constant. From this formula, one 
readily deduces that all moments of the energy (and therefore all moments of the amplitude of the oscillator) are of 
the order A which itself grows linearly with e, the distance from threshold. 

Another remarkable feature of the order parameter near threshold is that the time series t ^-^xit) exhibits on-off 
intermittency. Again this intermittent behaviour, first discovered in coupled dynamical systems [46| and in a system 
of reaction-diffusion equations [i^l, is beheved to be generic (see e.g., [11]) when an unstable system is coupled to a 
system that evolves in an unpredictable manner (multiplicative noise). 

Multiscaling and intermittency lead however to the following very puzzling problems : 

Problem 1 : The multiscaling behaviour given in equation (|18p is very clearly seen in computer experiments but 
does not seem to be observed in 'real' experiments. Similarly, It is surprising that, despite the genericity of the 
on-off intermittency mechanism (which is well established mathematically) this effect has scarcely been reported in 
experiments. One might expect that any careful experimental investigation of an instability should reveal on-off 



intermittency when the system is close to the onset of instability, and is hence sensitive to unavoidable experimental 
noise in the control parameters. 

Problem 2 : In well known works [13, HH, it was predicted, using perturbation theory, that the scaling of the 
stochastic bifurcation should be the same as that of the deterministic Hopf bifurcation, equation (|17p . This problem 
is analyzed in section HVl 




FIG. 9: Behaviour of the moments near the reentrant bifurcation threshold. The noisy double-well oscillator subject to 
parametric white noise is simulated with a fixed parameter a — 0.2. The critical value of the noise amplitude corresponding to 
the reentrant transition is found to be Ac — 0.627 : for a noise amplitude A belonging to the interval (0.627-1.70) the origin is 
stabilized by the noise. Even-order moments (x^), (a:*) and (x^) are plotted versus the distance to threshold e = (Ac — A)/ Ac 
(symbols). Dashed lines respecting a linear behaviour (a;^") oc e are drawn to guide the eye. Inset: the mean square position 
(measured in the stationary regime) is non-zero for A < 0.627 and A > 1.70, i.e., the stationary state is extended. 

Problem 1 is solved in [H, [l^. In these works, it is shown that although the noise strength controls the transition 
between the absorbing and the oscillatory states, the amplitude of the noise is not the only relevant parameter 
responsible for multiscaling and intermittency. Rather, these effects are due to the zero-frequency mode of the noise; 
in other words these effects depend on the Power Spectrum Density (PSD) of the noise and not on its overall amplitude. 
In order to identify which part of the PSD of the random forcing really affects the dynamics, one needs a random 
perturbation with a spectral density more complex than that of white noise or of Ornstein-Uhlenbeck noise. A useful 
type of noise [49] is the harmonic noise C(t) whose autocorrelation function is given by 

(C(i) C{t + r))s = Aexp(-27r77|r|) (cos(27rl7r) + ^ sin(27rf7|T|)) , (19) 

where A is the noise amplitude; the corresponding PSD is 

The value of the PDS at zero frequency is therefore given by 

S Aij/ [^(77^ + n')] . (21) 



For harmonic noise, the amplitude A and the zero frequency amplitude S can be tuned independently. 



In figure [TUl we plot the time series of a Duffing oscillator with a random frequency perturbation given by a harmonic 
noise C(^)- We observe that the zero frequency amplitude S is the pertinent parameter controlling the intermittent 
regime : on-off intermittency disappears when the value of S is lowered. For the simple first-order model, 

X = {a + C{t))x ~ , (22) 

it is possible to determine analytically the transition line between intermittent behaviour and non-intermittent be- 
haviour in the (5*, a) plane, using the cumulant expansion introduced by Van Kampen Intermittency occurs 
when 

< I < 1 . (23) 

This criterion is checked numerically in figure [TTl 

Finally, it is possible to prove that multiscaling is a consequence of the intermittent behaviour of the order parameter 
x{t) (19] , Therefore as the zero frequency noise amplitude is reduced, both on-ofF intermittency and multiscaling are 
suppressed and normal scaling (|17|) is recovered. 

This analysis solves Problem 1 by explaining why many experimental investigations on the effect of a multiplicative 
noise on an instability do not display on-ofF intermittency. If the noise is high-pass filtered, as often required for 
experimental reasons, then the regime of intermittent behavior disappears. This is the case for instance in [sol ]: a 
ferrofluidic layer undergoes the Rosensweig instability and peaks appear at the surface. The layer is then subject to a 
multiplicative noise through random vertical shaking. Close to the deterministic onset, the unstable mode submitted 
to a colored noise does not display intermittency. Another experimental obstacle to observe intermittency is the 
presence of additive noise that destroys the symmetry of the system. 
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FIG. 10: Temporal traces of the amplitude of a Duffing oscillator subject to harmonic noise of constant amplitude A — 0.05. As 
the value S of the zero frequency amplitude of the noise is decreased from top to bottom, intermittent behaviour is suppressed. 



IV. REMARKS ON THE PERTURBATIVE ANALYSIS OF A NOISY HOPE BIFURCATION 

Problem 2 questions the relevance of perturbative expansions for studying noise induced bifurcations. Consider, 
for example, the parametrically driven damped anharmonic oscillator that naturally appears in the study of many 
instabilities [s^. Such a system is described by the following equation : 

mx + mjx ^ {e + A^{t)) X - x^ , (24) 

where e is the control parameter and the modulation ^(t) is of arbitrary dynamics and statistics: it can be a periodic 
function or a random noise. For small driving amplitudes A, Liicke and Schank [20l | have performed a Poincare- 
Lindstedt expansion called the Poincare-Lindstedt expansion [H, IH, [5l| (see (sl] for applications to field-theory). 
They have obtained an expression for the threshold ec(A) (at first order in A). Their result has been verified both 
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FIG. 11: Intermittency phase diagram for the solution of (|22l) with harmonic noise. (A) : intermittent behavior, (o) : non- 
intermittent behavior. The full line is the transition curve predicted by (|23p. 

numerically and experimentally and is also in agreement with the exact result obtained for the Gaussian white noise 
(in this case a closed formula is available for ec(A) for arbitrary values of A). Another result obtained in [13, HH is 
the scaling of the moments near the threshold, 

= s„ (e - e,(A))" + O ((e - e,)"+i) , (25) 

where the constant s„ depends on ^(t) and on A. The moments have a normal scaling behaviour : (a;^") scales as (a;^)". 
The bifurcation scaling exponent is equal to 1/2 and is the same as that of a deterministic Hopf bifurcation. However, 
this expression does not agree with the results for random iterated maps, for the random parametric oscillator and 
with recent studies on On-Off intermittency [Tsl . [Tsl . [il ] . These works predict that the variable x is intermittent and 
that the moments of x exhibit anomalous scaling, 

(a;2") ~K„(e-ec) for all n>0, (26) 

i.e., all the moments grow linearly with the distance from threshold. This multiscaling behaviour, confirmed by 
numerical simulations for a Gaussian white noise, was derived using effective Fokker-Planck equations. The origin 
of the contradiction between equations and (PS)) lies in the divergences that appear in the Poincare-Linsdtedt 
expansion. This fact, identified in [2ll|, implies that the results of are valid only for noises that do not have low 
frequencies. 

A. A simple model 

The importance of low frequencies in the noise spectrum can be seen analytically in a a model technically simpler 
than equation (I24p . To simplify our discussion we work on the first order stochastic equation : 

x = {e + A^{t))x-x^, (27) 

that we already encountered in the previous section. Here, the noise ^(t) is a Gaussian stationary random process 
with zero mean value and with a correlation function given by 

mm)^m~t'\). (28) 

The power spectrum of the noise, which is non- negative thanks to the Wiener-Khinchin theorem jTi], is the Fourier 
transform of the correlation function 

/+00 r + co 

dtexp{iLut){^it)^{0)) = dtexp{iujt)V{\t\). (29) 
-OO J — OO 



Applying elementary dimensional analysis to equation (|27p . we obtain the following scaling relations: 

a;^ti/(2p)^ ^^t-^/\ e^A^^t-\ (30) 

The dimension of the noise ^ is so chosen as to render the power spectrum 'D^lu) a dimcnsionless function. 

We first remark that the presence of noise does not modify the bifurcation threshold; indeed the Lyapunov exponent 
of the system is given by A = e and therefore the bifurcation always occurs at ec(A) = 0. 

When ^(t) is a Gaussian white noise, the stationary solution of the Fokker-Planck equation corresponding to 
equation ([27|) is given by 



2 

X \ . , e 



where F represents the Euler Gamma- function. The bifurcation threshold is given by e = 0; for e < 0, the solution (j3ip 
is not normalizable : the stationary PDF is the Dirac distribution 5{x) localized at the absorbing fixed point a; = 0. 
For e > 0, the solution (PT|) is normalizable and the moments of x are given by (a;^") = ^Y^iy^(A^)" . In the vicinity 
of the threshold, e is small and we recover multiscaling : the moments scale linearly with e, i.e.. 



(X 



e(A2)"-iF(n). (32) 



B. Exact solution for arbitrary noise 

For an arbitrary noise ^(i), it is still possible to obtain an exact solution of equation ([27]) for all times. First, we 
shall derive some remarkable identities satisfied by the exact solution. 
We define the differential operator L as follows, 

£ = ^ - AC(t) . (33) 
Let us call y\{t) the solution of the adjoint equation £^t/i = 0, which is given by 



yi(t) =exp (^-A^ 



(34) 



Multiplying equation (|77|) by the function yi and taking average values, we obtain 

(yiZx) = {eyix - yix^) . (35) 

Integrating the left hand side of this equation by parts and taking into account the fact that yi is in the kernel of the 
adjoint operator C'' , we derive the following relation 

e=- - + — Inyix). (36) 

{yix) dt 

Similarly, we have 

e=(^') + ^(ln(2/ix)). (37) 

Dividing both sides of equation (P7)) by and using the auxiliary variable U — 1/x'^, we observe that equation P7)) 
becomes a linear first order stochastic in U : 

^U = l-{e + AOU. (38) 

This equation can be solved exactly for all time by using the method of variation of constants. Introducing the initial 
value a:(0) = A that has the dimensions 



A ~ ^ (39) 



we obtain the following explicit formula for x : 

\exp{et + ABt) 



x{t) = '7"^ ' " . (40) 

1 + 2A2 exp {2eu + 2AB„) du 



We have defined here the auxiliary random variable Bt : 



Bt= f au)du. (41) 
Jo 



Because ^ is taken to be a Gaussian random process, Bt is also Gaussian. 

C. Importance of low frequencies in the noise spectrum 

In order to determine the scaling on the moments of x, we must evaluate expressions of the type (exp ABt). Because 
Bt is a Gaussian random variable, this quantity is given in terms of the variance of Bt : 

(B?) - f f\iiumv))dudv ^ f fvi\u-v\))dudv= pi^^ ^HA.. (42) 



The last integral is well defined at w = (the time t introduces an effective low frequency cut-off for uj ~ 
The behaviour of (Bf) for < ^ cx) depends on the behaviour of "Diuj) at w ^ 0. The following two cases must be 
distinguished : 

(i) The spectrum of the noise vanishes at low frequencies, i.e., 2?(0) — 0. Because 2?(a;) is an even function of u, 
we suppose that 2?(w) ^ tJ^ for cj ~ (we disgard non-analytic behaviour of the power spectrum at the origin. Such 
non-analyticity would correspond to long tails in the correlation function of the noise) . 

(ii) The power spectrum of the noise is finite at w = 0, i.e., 2'(0) > 0. 

In case (i), the long time limit of equation (j42|) is readily derived and we obtain (by using the Ricmann-Lebesgue 
lemma) 

(B?)^ p^Hdo; ^^^^ ^^^^ (43^ 



The variance of Bt has a finite limit at large times. 

In case (ii) , the integral on the right hand side of equation (|42|) diverges when t ^ oo and its leading behaviour is 



iB?) = -J_^ ^^v{-)du^ViO)t. (44) 
The variance of Bt grows linearly with time in the long time limit. 

D. Behaviour of the moments 

We have seen that the behaviour of the variance of Bt in the long time limit crucially depends on the on the 
behaviour of the noise spectrum at low frequencies. When I?(0) > the noise £,{t) acts dominantly as a white noise 
and the multiscaling behaviour p2p is recovered when t — > cx). 

When ViO) = 0, the random variable Bt has a finite variance even when t ^ oo. Using equation (|40p and keeping 
only the dominant terms in the long time limit, we obtain 

m exp jet + ABt ) 

x{t) ~ . =^ ■ (45) 

' 2 /p exp (2eu + 2AB„) du 



Therefore, we have 



(t) ~ ( 2 ^ exp {~2e{t - u) + 2A(B„ - Bt)) du ^ " (46) 



The fluctuations of the random variable being of order one, the integral on the r.h.s. of this equation is dominated by 
the contribution of the linear term — 2e(t — u). The contribution is maximal in the region u t. Therefore, we have 

~ , (47) 

where the constant C„ depends a priori on the statistical properties of the noise. This equation predicts that subject 
to a noise without low frequencies, the amplitude x{t) exhibits a normal scaling behaviour identical to that of the 
deterministic case. 



V. CONCLUSION 



The stochastic oscillator is an ideal model to study the effect of random perturbations on a nonlinear dynamical 
system. Various effects of noise can be demonstrated on this simple model : noise can shift bifurcation thresholds, 
can create new phases by destabilizing (resp. stabilizing) stable (resp. unstable) fixed points, can induce reentrant 
behaviour. The relevant parameter that determines the long time behaviour of the system is the Lyapunov exponent 
of the underlying linearized dynamics. Exact results can be derived for white noise and for dichotomous Poisson noise. 
When the system is coupled to the smooth Ornstein-Uhlenbeck process, analytical treatments have to rely on various 
approximations, however qualitative aspects on time-correlations are well understood. 

The effect of noise on the scaling behaviour is subtle : intermittency and multiscaling can appear if the relative 
weight of the zero-frequency mode of noise is large enough. If low frequencies are filtered out, the noise-induced 
bifurcation becomes qualitatively similar to the deterministic transition. These features are difficult to extract from 
a perturbative analysis of the Langevin equation : they appear as divergences in the perturbative expansion which 
must be resummed in order to get correct results. Indeed, multiscaling of the order parameter can not be revealed 
from a finite order truncation of the expansion. 

Despite the formal simplicity of the problem, exact results have been derived only recently and many questions still 
remain to be addressed. For example, little is known about the relevance of the noise spectrum on barrier crossing 
problems and on stochastic resonance. Analytical results for higher dimensional systems and for stochastic fields are 
rare : such results could be relevant for the study of the trapping of quasi one-dimensional Bose-Einstein condensates 
in random potentials |53| • 
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